Persistent spin current in mesoscopic ferrimagnetic spin ring 
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Using a semiclassical approach, we study the persistent magnetization current of a mesoscopic 
ferrimagnetic ring in a nonuniform magnetic field. At zero temperature, there exists persistent spin 
current because of the quantum fluctuation of magnons, similar to the case of an antiferromagnetic 
spin ring. At low temperature, the current shows activation behavior because of the field-induced 
gap. At higher temperature, the magnitude of the spin current is proportional to temperature T, 
similar to the reported result of a ferromagnetic spin ring. 
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INTRODUCTION 



Persistent charge current in a mesoscopic metal ring 
was predictedi and observed^ a decade ago. In such a 
ring threaded by a magnetic flux, if the phase coherence 
length of electrons is larger than the size of the ring, 
then the electrons can pick up an Aharonov-Bohm (AB) 
phase after circling the ring once. Such a phase lag (or 
advance) would lead to a persistent current, which is a 
periodic function of the threaded magnetic flux,'^ and can 
be detected via the magnetic response of the (isolated) 
ring. The phase lag (or advance) can also be of geomet- 
ric origin (Berry phase) 4 It has been proposed that a 
Berry phase can appear for an electron moving around 
the metal ring that subject to a textured magnetic field 
(or magnetization).^ This geometric phase, which de- 
pends upon the solid angle associated with the textured 
magnetic field, can lead to persistent charge and spin 
currentsi^ A similar geometric phase appears due to the 
spin-orbit interaction in one-dimensional ringSjSii which 
is a manifestation of the Aharonov-Casher (AC) effect.^ 
More studies on the persist ent cu rrent related to the AC 
effect can be found in Refs. lOllollTll 

With the advance of spintronics and quantum 
computation,^^ it becomes more important to under- 
stand the behavior of the spin current. Among these in- 
vestigations, spin transport in pure spin systems plays a 
special role since there is no complication from charge de- 
grees of freedom. In a recent paper, using a semiclassical 
spin wave analysis, Schiitz et al. predicted the existence 
of persistent spin current in a mesoscopic ferromagnetic 
(FM) spin ring in a nonuniform magnetic fieldji^ The 
FM spin ring being considered is a charge insulator with 
Heisenberg spin interaction, and the spin current is car- 
ried by magnon excitations. Similar to the case of charge 
transport in a metal ring subject to a textured magnetic 
field)^ the magnon in a mesoscopic FM spin ring acquires 
a geometric phase from the (nonuniform) spin texture of 
the classical ground state. The persistent current is found 
to be zero at temperature T = 0, and proportional to T 
when fcsT is larger than the field- induced energy gap of 
the magnons. 



Similar method has been applied to an antiferromag- 
netic (AFM) spin ring with a Haldane gap.^^ As com- 
pared with the FM case, there are some subtleties in 
using the semiclassical method in the AFM case. Due 
to the problem of infrared-diverging magnetization, the 
spin-wave approach is not valid for AFM spin chains 
with half-integer spins.^^ That is the reason why only the 
integer-spin cases are considered in Ref. |3 Nonetheless, 
in the integer-spin case, an additional staggered field in 
the direction of the classical magnetization vectors still 
has to be introduced. Its value needs to be determined 
self-consistently before quantitative predictions can be 
made. The authors of Ref. il4l find that, unlike the case 
of the FM spin ring, the persistent spin current in an 
AFM spin ring can be nonzero at T = due to quantum 
fluctuations. When the spin correlation length is much 
longer than the size of the ring, the magnitude of the spin 
current exhibits sawtooth variation with respect to the 
geometric phase, similar to the case of persistent charge 
current in a metal ring. Recently, the investigation has 
been extended to an anisotropic FM spin ringfi^ a spin- 
1/2 AFM spin ringfiSiii and an anisotropic AFM spin 
ring.^^ 

In this paper, we study the persistent spin current in 
a ferrimagnetic (FIM) spin ring with alternating spins 
S'^ and under a textured magnetic field. Contrary 
to the AFM case, the problem of infrared-diverging mag- 
netization does not exist in the present FIM case, no 
matter whether the constituent spins are integer or half- 
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Thus the self-consistently determined 



staggered field needs not be introduced, and physical 
quantities can be calculated directly as long as system 
parameters are known. We find that the FIM spin ring 
can have either FM or AFM characteristics. For exam- 
ple, a quantity proportional to — plays a role 
similar to the Haldane gap in the AFM spin ring. More- 
over, a nonzero spin current exists at T = 0, again similar 
to the case of the AFM spin ring»i4 On the other hand, 
when the thermal energy is higher than the field-induced 
gap, the magnitude of the spin current is proportional to 
temperature T, similar to the case in the FM spin ring. — 
This paper is organized as follows: We present the spin- 
wave analysis in Sec. II. The results of numerical calcula- 



2 



tions are shown in Sec. Ill, and Sec. IV is the conclusion. 



II. THEORETICAL ANALYSIS 



The Hamiltonian of the ferrimagnetic Heisenberg spin 
ring in a nommiform magnetic field hj = g^sB^fj) is 



H = J 



E 

jeAuB 



- E 

]ieA,j2'EB 



■ S: 



(1) 

where J > 0, and the index j refers to one of the alter- 
nating ji, j2 sites. That is, sublattice-A can be labelled 
either by j (an odd integer) or ji = (j + l)/2; sublattice- 
B can be labelled by j (an even integer) or j2 = j /2. 
There are even number of lattice sites N . The length of 
the ring is L and the lattice spacing a = L/N . Periodic 
boundary condition Sj+N = Sj is used. On the classical 
level, Sj are replaced by classical vectors Srhj . The clas- 
sical ground state {rhj\ can be determined from angular 
variations with respect to each rhj , which give 



(2) 



where Xj are Lagrange multipliers. It shows that the 
magnetization aligns parallel to the sum of external and 
exchange field, as expected. If the Zeeman energy is much 
smaller than the exchange energy between spins, then 
and w^j^ would be nearly antiparallel to each other, as 
shown in Fig. 1(a). Moreover, due to nonzero magneti- 
zation in the present FIM model, would lie nearly 

along the direction of hj , instead of nearly perpendicular 

to hj as in the AFM casein (we take S"^ > in this 
paper). 

When quantum fluctuations are considered, following 
the treatment introduced by Schiitz et al.^^^^ each spin 
operator is expanded using unit vectors that form an or- 
thogonal triad {ej 



(3) 



where = ej±ie|. We focus on systems with large spins 
(which could be integer or half- integer), and introduce 
the Holstein-Primakoff bosons aj , hj , 
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Substituting these into Eq. (1), we obtain H 
ff-L + H' with 
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FIG. 1: (a) Classical spin configuration of a FIM spin ring 
in a crown-shaped magnetic field, (b) Relative orientation 
for adjacent local triads. The condition e? — e?ii = rh, x 
rhj+i/\mj 



The condition e_, — Sj^i 
X rhj+il is imposed. The circle in the middle is the 



equator of the sphere. 



= E {sf,^-[jS''{m,,^i+ rh,, ) - h 



neA,j2eB 



+S. 



'Bl^ 



[jS'^irhj, -l-mji+i) - hj2 I 



-O(V^). 



(7) 



i?" and are of the orders of 0{S^) and 0{S) re- 
spectively. To the order of 0{S), H' is zero because 
of Eq. (2) and will be neglected in the following. In 
order to simplify the Hamiltonian, we choose the local 
triads with the connection shown in Fig. 1(b), in which 
= e|_|_]^ = rhj X mj+i/\rhj x rhj+i\ (this is referred 
to as a choice of gauge). Such parallel transported tri- 
ads are related to the original local triads by local rota- 
tions (gauge transformations), — e^^'^^^^+^ e^, where 
ujj^j^i is the angle of rotation around rhj that takes e| 



to e 



£2 23 



Using the new triads in the Hamiltonian, we obtain 
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J2 i^'n-n-n + hf/,,b,X 



(8) 



Jl J2 



where Hr 



j '"3 



T.,,,,AS''hf^+S^hfJ with hf^ EE h,, ■ fn,, and 
hf^ EE hj^ -rhj^. Also, 



31^2 



+ /i.e.} . (9) 



As long as the ring is not too small, neighboring spin vec- 
tors would be nearly antiparallel to each other [as shown 
in Fig. 1(a)] withmj-TOj+i = -l + 0(l/iV). This further 
simplifies the Hamiltonian to be 



31,32 



(10) 



The hopping of boson-a from site-ji to site-(ji + 1) ac- 



ji+l^i2y 



quires aphase (wji^j2+'^j2^ji)~(^j2^ji+i 
After circling around the ring once, the boson gains a cu- 
mulative phase = X^jLi ("l)"'"''^ ('^j^j+i 
which is the holonomy angle of the parallel transport 
and equals the solid angle extended by the classical spin 
texture \mj\?^ It can be shown that 
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^ j2 ^32 



-31+^) 



(11) 
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which is a gauge-invariant expression because the e 
vectors always appear in pairs (note that e^^2+i = sT)- 
Therefore, this cumulative geometric phase is indepen- 
dent of the choices of the local triads. 

The phase factors in Eq. (10) can be made implicit 
by merging with the boson operators. Such new boson 
operators would satisfy the twisted boundary conditions: 

With the help of 



the transformations, 



a-k 



/|^e---a,, &.^/|E^""^'^.' (12) 

3&A jeB 



in which 

2iT f n\ N 

kr, = — n H , n = 0, 1, 2, ... 1 , 

L V 27r y ' ' ' ' 2 

to conform with the twisted boundary condition, the 
Hamiltonian becomes 

H = Hr 



- E {2JS^ + h^) 4+,a, + (2 + /if) 

-2 JVS-^^s E (afe&fc + fltfotA cos(/ca), (13) 



where we have assumed that the applied magnetic field 
(as well as {rhj}) has only one Fourier component with 
momentum q to simplify the expression. A crown-shaped 
magnetic field with azimuthal symmetry has the g = 
component only. For convenience, we consider the crown- 
shaped magnetic field below. For a large FIM spin ring 
in a weak magnetic field, we also have = —h^ = —ho 
[see Fig. 1(a)] with ho being positive . With the help of 
the Bogoliubov transformation. 



= ak cosh 9k + PI sinh 0k ; 
= ak sinh 9k + pl cosh 9k , 



and choosing 



tanh(26'fc) 



2\/SAS^cos(A:a) 



the Hamiltonian is finally diagonalized as 



(14) 



(15) 



H 



E 



k I + o I + 



NJS^ 



■(1 + 7), 



(16) 



where the two energy branches are 



y^(l-7)2 4-47sin2(/ea)± (1 - 7) 



(17) 



with 7 = / < 1. As mentioned before, contrast to 
the AFM case^i^ the staggered field needs not be intro- 
duced in the present FIM case. Thus physical quantities 
can be calculated directly as long as the system parame- 
ters are known. 

Similar to the case of a FIM spin chain under a uni- 
form magnetic field, the magnons with energy (e^) 
in the present case correspond to the ferromagnetic (anti- 
ferromagnetic) excitations. The energy gaps of these two 
branches are = and Cq = 2JS^[1 — 7) — /ig, re- 
spectively. That is, a gap is induced by the applied field 
for the ferromagnetic excitations, while the gap of the 
antiferromagnetic excitations is reduced by the applied 
field. In the absence of external magnetic field, the fer- 
romagnetic branch ejT becomes gapless with quadratic k- 
dispersion at small k, which corresponds to the Goldstone 
mode due to the spontaneously broken rotational symme- 
try. Calculations using quantum Monte Carlo method 
yields nearly the same curve for e^, but is separated 
from €^ with a larger (fc-independent) gap^S*^ Such a 
discrepancy is reduced when the spins are larger and the 
semiclassical formalism works better 

Since the magnons are non-interacting, it is straight- 
forward to obtain the free energy (ks = 1) 



k 



4sinh ( ^ 



(I+7). 



sinh 



(18) 



Furthermore, it can be explicitly shown that the longitu- 
dinal (gauge-invariant) spin current 



dF{n) 



(19) 



similar to the relation for persistent charge current in a 
normal metal ring, but with f2 replacing magnetic flux 
j^ji3ji4 Pj-qj^ ^j^jg relation we obtain the magnetization 
current. 



where 



g^J■B 



Idel 
h dk 



9MB 
L 



fc, a— ± 



2JS^ 



7sin(2fca) 



\J{1- lY + 47sin^(fca) 



(20) 



(21) 



are the velocities of the magnons, and — 
1/ [exp (ej?/r) — 1] are the Bose occupation numbers. 



III. BEHAVIOR OF THE PERSISTENT 
MAGNETIZATION CURRENT 

For the FIM case, the magnetization current at T = 
is nonzero even for vanishing magnon numbers (see 
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FIG. 2: Amplitude of the magnetization current at zero tem- 
perature as a function of 7, plotted for three different ring 
sizes. 



Eq. 20), similar to the AFM case, 

Q 7sin(2fca) 



m / J 



\J{1- lY + 47sin2(fca) 



(22) 



where = -{2gfiB /fi){JS'^ /N). The magnon veloc- 
ity within the summation is a periodic function of k. 
Therefore, after summing over the first Brillouin zone, 
the current is zero if the fc-points are distributed sym- 
metrically (fi — tt). For other values of fl, the sum- 
mation is nonzero but the magnitude of the magnetiza- 
tion current decreases rapidly as 1 — 7 becomes larger. 
Compa ring with th e AFM case^^ it can be seen that 
A = y^{l — 7)^/47 plays a role similar to the Haldane 
gap, and its inverse determines the scale of the spin cor- 
relation length ^. Therefore, when 7 is small enough 
such that A ^ 2tt/N, the spin correlation length <^ L; 
while if 7 ~ 1, such that A <C 2tt/N, we have ^ ^ L. 
Therefore, by varying the ratio of the two different spins, 
qualitatively different regimes, L and ^ ^ i, can be 
reached. In Fig. 2, the maximum amplitude of the mag- 
netization current /max is plotted as a function of 7. The 
functional form of /max (7) shows a very clear crossover 
between these two different regimes. 

Thermal energy could excite magnons and generate 
larger magnetization current. Typical influence of the 
temperature on the magnitude of the magnetization cur- 
rent can be seen in the inset of Fig. 3. We have also 
studied the dependence of /max on the parameters T, ho, 
and 7. From the magnon dispersion relations in Eq. (17), 
we expect activation behavior for /max(T') at low temper- 
ature T < Cq = ho. When the thermal energy is larger 
than the fleld- induced energy gap ho, /max(T') should be 
proportional to T, similar to the behavior of the persis- 
tent spin current in a FM spin ring.^-i These behaviors 
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FIG. 3: Amplitude of the magnetization current as a function 
of temperature T. Solid (dotted) line corresponds to a system 
(iV = 100, 7 = 0.9) in a magnetic field ho/2JS^ = 0.01 (0.05). 
Other curves with different 7 aU have ho/2JS^ = 0.01. Inset: 
The variation of the magnetization current with respect to the 
change of Q. and T (N — 100, 7 — 0.8). The temperatures for 
the curves with the smallest amplitude (indiscernible from 
the horizontal axis) to the largest amplitude are T/JS^ = 0, 
0.005, 0.010, 0.015, 0.020, and 0.025. 



IV. CONCLUSION 

We extend the work of Schiitz et ali^^ and stud- 
ied the persistent spin current in a FIM spin ring. At 
T = 0, the functional form of the magnetization current 
Im{l) shows distinctive behaviors above/below a thresh- 
old value of 7, which depends on the size of the ring. 
When thermal excitation can overcome the field-induced 
energy gap = ^Oi the magnitude of the spin current 
grows linearly with temperature T. 

Both the persistent charge current in a metal ring and 
the persistent spin current in a spin ring are related to ge- 
ometric phases (albeit of different origins). As in the case 
of the FM spin ring, the induced electric voltage would 
be of the order of nV, which poses a stringent experimen- 
tal challenge. Recently, several types of FIM spin-chain 
compounds have been synthesizediSL2& To obtain a FIM 
spin chain with 7 > 0.8, one might need to introduce 
rare earth elementsj^S or to fabricate a ring composed of 
magnetic molecules with large spins 

On the theoretical side, several questions remain open, 
such as the generalizations to quasi-lD spin rings or spin 
rings with itenerant electrons. The realistic effects of 
disorder, interaction, or contact leads remain unknown 
in the spin ring case. Further studies are desired. 
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can be seen in Fig. 3. For the upper two curves with 
7 = 0.9, there is a significant amount of spin current at 
zero temperature and the activation behavior is implicit. 
Also, the persistent magnetization current at very low T 
is independent of the value of ho (see Eq. 22). 
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